Conversion of terahertz wave polarization at the boundary of a layered 
superconductor due to the resonance excitation of oblique surface waves 
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We predict a complete TM-f-^-TE transformation of the polarization of terahertz electromagnetic 
waves reflected from a strongly anisotropic boundary of a layered superconductor. We consider the 
case when the wave is incident on the superconductor from a dielectric prism separated from the 
sample by a thin vacuum gap. The physical origin of the predicted phenomenon is similar to the 
Wood anomalies known in optics, and is related to the resonance excitation of the oblique surface 
waves. We also discuss the dispersion relation for these waves, propagating along the boundary of 
the superconductor at some angle with respect to the anisotropy axis, as well as their excitation by 
the attenuated-total-reflection method. 
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Surface waves play a very important role in many fun- 
damental resonant phenomena in different fields of mod- 
ern physics. A well-known example of such phenomena 
arc the Wood anomalies in the refiectivity and transmis- 
sivity of periodically-corrugated metal and semiconduct- 
ing samples (see, e.g., Refs. [iH^)- A more recent ex- 
ample is the extraordinary transmission of light through 
metal films perforated by subwavelength holes The 
excitation of surface plasmons can also result in an "in- 
verse" effect of resonant suppression of light transmission 
through perforated metal films with thicknesses less than 
the skin-depth 0. 

Recent interest on the effects listed above is also partly 
due to their possible applications, e.g., in photovoltaics, 
as well as the detection and filtering of radiation in the 
far-infrared and visible frequency ranges. From this point 
of view, it is interesting to consider layered superconduc- 
tors, instead of metals, for designing devices which op- 
erate using the excitation of surface waves. Indeed, the 
characteristic frequencies of surface waves in these ma- 
terials belong to the terahertz frequency range, which is 
important for different applications, but is still hard to 
reach for electronic and optical devices. 

Layered superconductors are either artificially-grown 
stacks of Josephson junctions, e.g., Nb/Al-A10a;/Nb, 
or natural high-temperature superconductors, such as 
Bi2Sr2CaCu208-i-5. These materials contain thin super- 
conducting layers separated by thicker dielectrics. Exper- 
iments on the c-axis transport in layered superconductors 
justify the use of a theoretical model in which the super- 
conducting layers are coupled by the intrinsic Josephson 
effect through the layers (sec, e.g., Refs. [^. Due to 
their layered structure, these superconductors possess a 
strongly-anisotropic current capability. 

The multi-layered structure of Bi2Sr2CaCu208+d' (and 
similar superconductors) supports the propagation of 
specific Josephson plasma electromagnetic waves ( JPWs) 



(see, e.g., the review [10| and references therein). The 
spectrum of JPWs lies above the so-called Josephson 
plasma frequency, wj ^ 1 THz. In a semi- infinite sample, 
apart from bulk JPWs, surface Josephson plasma waves 
(SJPWs) can propagate along the interface between an 
external dielectric and a layered superconductor. These 
waves are similar to the surface plasmon-polaritons in 
normal metals. As shown in Refs. I 111 , the spectrum of 
SJPWs also lies in the terahertz range and consists of 
two branches: one above ojj and another below it. Simi- 
larly to the optical resonance phenomena in normal met- 
als, the resonance excitation of SJPWs is accompanied 
by Wood anomalies in the refiection and transmission of 
terahertz waves through layered superconductors. How- 
ever, the strong anisotropy of the current capability can 
result in new resonance phenomena, specific to layered 
superconductors. 

In this paper, we predict one of such unusual resonance 
phenomena. We show that the resonance excitation of 
surface waves can be accompanied by a complete con- 
version of the polarization of the terahertz radiation af- 
ter reflection from an anisotropic boundary of a layered 
superconductor. For example, if the incident wave has 
a transverse magnetic (TM) polarization (a wave with 
the magnetic field parallel to the sample surface), the re- 
flected wave is completely transverse electric (TE) (with 
the electric field parallel to the sample surface), at an 
appropriate choice of the direction of the incident wave 
vector. We emphasize that this phenomenon is not an 
analog of the Brewster effect (where light with TM polar- 
ization, for a definite incidence angle, is perfectly trans- 
mitted through a transparent dielectric surface, with no 
reflection, and, therefore, the reflected light only con- 
tains the TE polarization). However, in the phenomenon 
predicted here, the reflected wave has a TE polarization 
when the incident radiation does not contain the wave 
with this polarization, i.e., when the incident wave is com- 
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Figure 1: (Color online) Schematic geometry of the problem 
for the excitation of oblique surface waves. Here k', k*^, and 
K are the wave vectors of the incident, reflected, and oblique 
surface waves. Note that S=superconductor and I=insulator. 



pletely TM polarized. In other words, the reflected wave 
is TE-polarized due to the conversion of the incident TM 
radiation, but not due to separation of the TE wave from 
the incident mixed (TE+TM) radiation. 

This polarization conversion has a clear physical expla- 
nation. We study here the so-called "attenuated-total- 
rcflection" (ATR) method for producing surface waves. 
The surface wave is excited by the polarized radiation 
(with, e.g., TM polarization) incident on the supercon- 
ductor from a dielectric semi-spherical prism separated 
from the sample by a thin vacuum gap (see Fig. [T]). 
We assume that the crystallographic c-axis of the super- 
conductor lies in the sample surface, and, therefore, the 
boundary of the superconductor is strongly anisotropic. 
For arbitrary direction of the incident radiation, the wave 
vector K of the excited surface wave is oriented at some 
angle d with respect to the anisotropy axis c, and we 
call these waves oblique surface waves (OSWs) [121 ■ Due 
to the anisotropy, the OSWs contain all components of 
the electric and magnetic fields and, therefore, the radi- 
ation reflected from the superconductor contains waves 
of both, TM and TE, polarizations. In other words, due 
to the anisotropy, the reflected wave is, generally, un- 
polarized. Note that the ratio between the amplitudes 
of the reflected TM and TE waves is controlled by the 
parameters of the problem, such as the wave frequency 
w, angles -d and (p, the thickness h of the vacuum gap, 
etc. The most important issue is the possibility to choose 
the angles ■(? and if for the incident radiation in such a 
way that the amplitude of the TM reflected wave van- 
ishes. Recall that, in the isotropic case, one can choose 
the optimal value of the gap thickness h to provide the 
complete suppression of the reflected wave under reso- 
nance conditions. Similarly, in the anisotropic case, we 
can choose the optimal value of the angle to provide 
the complete suppression of the reflected TM wave under 
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Figure 2: (Color online) Geometry for oblique surface waves. 



resonance conditions. In this case, the desired complete 
conversion (TM — > TE) of the polarization takes place. 

Below wc discuss the dispersion relation for oblique 
surface waves, obtain their excitations by the ATR 
method, and calculate the reflection coefficient Rtm for 
the TM wave and the conversion coefficient Ttm-^te for 
the mode conversion from TM to TE. Similar conver- 
sion phenomena can be observed for transitions from TE 
to TM modes, from incident ordinary waves to extraor- 
dinary ones, and from extraordinary waves to ordinary 
ones. 

Oblique surface waves. — Consider a semi-infinite lay- 
ered superconductor with the crystallographic c-axis par- 
allel to the sample surface (see Fig. [5]). We choose the 
coordinate system in such a way that the z-axis coincides 
with the crystallographic c-axis of the superconductor 
and the y-axis is perpendicular to the sample surface. 
We study oblique surface waves, propagating along the 
interface vacuum-layered superconductor at an angle 
with respect to the c-axis, with electric E and magnetic 
H fields in the form, 

E(a:, y, z, t), H(a;, y, z, t) oc exp[i(Kp — ujt)], 

where p = (x, z) is the radius vector in the xz-plane, 
— {kx,kz) is the wave vector of the oblique surface 
wave, |k| > w/c, and c is the speed of light. It is suit- 
able to present the electromagnetic field as a sum of two 
terms that correspond to the ordinary and extraordinary 
evanescent waves, 

E(2/) = E(°)(y)+E(°)(y), H(y) = H(°)(y)+H(°)(y). (1) 

The electric field of the ordinary wave and the mag- 
netic field of the extraordinary wave are orthogonal to 
the anisotropy z-axis: 



E(°)(y)=H(°)(y)=0. 



(2) 



For each of these waves, the Maxwell equations give the 
relations between the components of the electric E™'^ and 
magnetic H™'^ fields in the vacuum and the same expo- 
nential law for the decay of both waves along the y-axis. 



E(°i(y), EtUy), H(°i(y), Hi^) (y) ex exp(-»fc-y) 



(3) 
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with imaginary normal component of the wave vector, 

k-^^{kl^kl-u?|c^fl\ (4) 

The electromagnetic field inside the layered supercon- 
ductor is determined by the distribution of the gauge- 
invariant phase difference x of the order parameter be- 
tween neighboring layers. This phase difference can be 
described by a set of coupled sine-Gordon equations (see 
review [lol | and references therein). In the continuum 
and linear approximation, x can be excluded from the 
set of equations for electromagnetic fields, and the elec- 
trodynamics of layered superconductors can be described 
in terms of an anisotropic frequency-dependent dielectric 
permittivity with components Sc and Soh across and along 
the layers, respectively jl3l |. 
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(5) 



Here we introduce the dimensionless parameters f2 = 
w/wj, Vah = 47r(Ta{,/ea;j7^, and = 47r(Tc/ewj, 7 = 
Ac/Aafo ^ 1 is the current-anisotropy parameter, Ac = 
c/Loje^/"^ and Aa& are the magnetic-field penetration 
depths along and across the layers, respectively. The 
relaxation frequencies Vab and Vc are proportional to 
the averaged quasi-particle conductivities Gab (along the 
layers) and CTc (across the layers), respectively; ujj = 
(STreDjc/heY^^ is the Josephson plasma frequency. The 
latter is determined by the critical Josephson current 
density jc, the interlayer dielectric constant e, and the 
spatial period D of the layered structure. 

Contrary to the vacuum, the Maxwell equations give 
different laws for the decays of the ordinary and extraor- 
dinary waves in layered superconductors, 

E(°)(y), H(°)(y)cxexp(zA:(°)y), 
E(^)(2/), H(^)(y)(xcxp(zfc(^)y) (6) 

with 

ki°^=i{kl + kl^cj'eab/cy^', (7) 

41' =^{kl + klejeab ~ w'sc/c")'/' ■ (8) 

From the continuity conditions for the tangential com- 
ponents of the electric and magnetic fields at the interface 
y = 0, we derive the dispersion equation of the OSWs, 
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= 0. (9) 




Figure 3: (Color online) Dispersion curves Q,{cn/ujj) of the 
OSWs for e = 16, 7 = 200, and Uab = Vc ^ Q. Curves 1-6 
correspond to i9 = 0°, 10°, 20°, 30°, 60 °, and 90 °, respec- 
tively. The hollow circles at f2 = Sli = e/(e — 1) show the 
end-points of the dispersion curves. In these points, the ex- 
traordinary wave, which is part of the OSW, is transformed 
into a bulk wave. 



Figure [3] illustrates the numerically calculated dispersion 
curves of the OSWs for different values of the angle d. 

Excitation of the OSWs and the polarization 
conversion. — Here we investigate the OSWs exci- 
tation using the ATR method. Consider a TM polarized 
wave incident from a dielectric prism with permittiv- 
ity Ep onto a layered superconductor separated from 
the prism by a vacuum gap of thickness h (see Fig. 
[1]). The incident angle ip exceeds the limit angle 

— 1/2 

9t ~ arcsin(£p ) for the total internal reflection. In 
this case the incident wave can excite a surface wave if 
the resonance condition, loe^J"^ sin(y5 = ck is satisfied. 

The electromagnetic field in the dielectric prism is a 
sum of three terms that correspond to the incident TM 
polarized wave and two reflected waves with the TM and 
TE polarizations. Thus, the x-component of the electric 
field in the prism can be presented as 

EP{x,y,z,t) = [Elcxpiikpyy) 

+ (-E^Jtm +-E'^TE)cxp(-ifcpj^2/)]exp[i(/«p- wt)] (10) 

where El, E^.^,^, and i?^xE amplitudes of the 

x-components of the electric field for the incident waves 
and for the reflected TM and TE waves, respectively; 
kpy = {oj/c) Ep^ ^ cos tp is the normal component of the 
wave vector of the incident and reflected waves. Other 
components of the electromagnetic field in the prism are 
expressed via the amplitudes E]. , E^ rpjyj, and E^. using 
the Maxwell equations. 

In the vacuum gap, we deflne the electromagnetic field 
as the superposition of four evanescent waves. They are 
two TM modes with imaginary y-components ±fc™'^ of 
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Conclusions. — We predict the complete polariza- 
tion conversion for terahertz waves reflecting from the 
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Figure 4: (Color online) The reflectivity -Rtm and the trans- 
formation coefficient Ttm^te, both versus ip. 



the wave vectors given by Eq. dU and two TE modes with 
the same wave vectors. In the layered superconductor we 
have, similarly to the previous section, ordinary and ex- 
traordinary waves with the normal components of the 
wave vectors given by Eqs. ([7]) and ([8]). Then, using the 
continuity boundary conditions for the tangential compo- 
nents of the electric and magnetic fields at the interfaces 
dielectric-vacuum and vacuum-layered superconductor, 
we obtain eight linear algebraic equations for eight un- 
known amplitudes (for 4 waves in the vacuum, 2 waves 
in the layered superconductor, and amplitudes i^^TMi 
^xTE of the reflected waves in the dielectric prism). 
Solving these equations, one can derive the reflectivity 



i\2 



coefficient i?TM(V;^) = |EtmI /|E 
for the TM wave and the coefficient Ttm- 
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of the trans- 



formation from the TM mode to the TE one. Figure H] 
shows the dependences of these coefficients on the angle 
ifi, calculated for Sp = 10, £ = 16, 7 = 200, n = 1.00035, 
i^ab = Vc = 10~^, h = c/2wj, i9 = 15.6°. The curves 
demonstrate the complete transformation of the TM- 
polarized wave to the TE wave, for the optimal value 
of the angle d = 15.6°, under resonance condition. The 
positions of the minimum in the Rtm (</') dependence and 
of the maximum in the Ttm->te('/') dependence coincide. 
This position corresponds to the most effective excitation 
of the oblique surface wave. Similar effects can be also 
observed for the TE— ;>TM complete transformation of the 
polarization, as well as for the complete transformation 
of the incident ordinary waves to extraordinary and vice 
versa. Figure [5] shows the dependences of the reflection 
coefficient i?TE('(f, "(9) and the coefficient TTE-i-TM('if , i?) 
of the transformation from the TE mode to the TM one 
calculated for = 10, e = 16, 7 200, Vl = 1.0004, 
Vah = v^ = 10-3, h ^ 0.18c/w,7, and ?9 22°. 
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Figure 5: (Color online) The reflectivity -Rte and the trans- 
formation coefficient rTE->TM, both versus (p. 



strongly-anisotropic surface of layered superconductors. 
The origin of this unusual effect is similar to the Wood 
anomalies known in optics, and is related to the reso- 
nance excitation of the oblique surface waves. 
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